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We derive an exact lower bound to a universal measure of frustration in degenerate ground states of quantum 
many-body systems. The bound results in the sum of two contributions: entanglement and classical correlations 
arising from local measurements. We show that average frustration properties are completely determined by the 
behavior of the maximally mixed ground state. We identify sufficient conditions for a quantum spin system to 
saturate the bound, and for models with twofold degeneracy we prove that average and local frustration coincide. 
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Many body systems are typically modeled by Hamiltoni- 
ans that are sums of local terms. Each local term operates 
only on a part of the entire system and acts to minimize the 
corresponding energy. If different subsystems overlap, the 
competition among the different local terms can preclude the 
existence of configurations satisfying all such minimizations 
simultaneously, a phenomenon known as frustration UHsl. 
For classical Hamiltonian systems, frustration is associated to 
some nontrivial geometric property of the system itself Q3J] . 
On the other hand, due to quantum non-commutativity and 
entanglement , classically unfrustrated systems may ad- 
mit frustrated quantum counterparts [4, 6-8J. The existence 
at a qualitative level of a relation between frustration, non- 
commutativity, and entanglement has motivated recent efforts 
aiming at qualifying and quantifying frustration in the quan- 
tum domain. Frustration criteria, reducing to the Toulouse 
conditions in the classical case |Q]], have been introduced re- 
cently for systems with non degenerate ground states fl, and 
a quantitative relation between a universal measure of frustra- 
tion and ground-state bipartite entanglement has been estab- 
lished in the form of an exact inequality. This result cannot 
be extended unambiguously to degenerate ground states, be- 
cause of the existence of infinitely many pure ground states 
with different entanglement properties. In addition, the pos- 
sibility for mixed ground states introduces a further element, 
classical statistical correlations, in the quantification of frus- 
tration at the quantum level. 

The present work introduces a theory of frustration for 
quantum systems with arbitrarily degenerate ground states. 
While including the one discussed in Ref. [9] as a special 
case, it provides a general picture according to which the in- 
terplay between degeneracy and superselection rules results 
in a highly nontrivial relation between frustration and differ- 
ent types of correlations, classical and quantum. It provides 
a unified treatment and a rigorous quantification of frustra- 
tion, expressed in terms of a universal inequality. The latter 
assesses the relative weights of geometry, bipartite entangle- 
ment, and shared classical correlations, identified as the three 
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fundamental sources of frustration in quantum and classical 
systems. The unified approach is based on the concept of 
maximally mixed ground state (MMGS), that is the statisti- 
cal mixture with equal a priori weights of all possible pure, 
degenerate ground states. This state always satisfies all the 
symmetries of the corresponding quantum Hamiltonian and 
therefore stands as the natural candidate to provide the rele- 
vant information on the global characteristics of a quantum 
system. Thus the concept of MMGS allows to introduce a 
complete classification of many-body quantum Hamiltonians, 
well beyond the limiting cases discussed in Ref. [9J]. Finally, 
we determine sufficient conditions for the saturation of the in- 
equality that constitute a further quantum generalization of the 
classical Toulouse criteria. 

Let us consider a many body Hamiltonian Hj = Y*s n s sum 
of local Hamiltonians lis- Frustration is the impossibility for 
the ground state po of the total Hamiltonian Hj to be entirely 
projected in the ground space of every local Hamiltonian h$, 
that is the local ground space. To be more specific, let us 
consider a total system T, associated to the Hamiltonian Hj, 
bipartite into a subsystem S , corresponding to the Hamilto- 
nian hs, and the remainder R. Define then p$ = tr«(pc), the 
reduced ground-state density matrix of S , obtained by tracing 
out all the degrees of freedom of R in pc- Finally, introduce 
n,5 , the projector onto the local ground space of subsystem S . 
A natural way to quantify frustration is to consider the dis- 
tance between pc and the projector onto the ground space of 
the local interactions JjJ: 

fs = 1 - tr(p G • Il s ® 1«) = 1 - tr s (psTls) . (1) 

This quantity measures how much the global ground state fails 
to accommodate for the ground states of the local interactions. 
Let us define d as the rank of n,s , namely the degeneracy of 
the local ground space associated to subsystem S . The Cauchy 
interlacing theorem II 1 Ofl yields that the universal measure of 
frustration fs is bounded from below by the first d largest 
eigenvalues of p s , arranged in descending order {Aj(p s )); 19]: 

f s >ef, 4« = l-^A|(p s ). (2) 

For non degenerate, pure ground states pc = l^cX^cU the 
quantity ei coincides with the distance of pc from the set 
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of states with Schmidt rank less or equal than d [11]. Such 
distance is not a full bipartite entanglement monotone: it van- 
ishes for all separable states and is non-increasing under local 
operations and classical communication (LOCC) [11]; how- 
ever, it vanishes also for all entangled states with Schmidt rank 
Dschmidt ^ d. It is then restored to a full entanglement mono- 
tone for all entangled states with Dschmidt > d. Therefore, 
for d — 1 the quantity ei is a faithful ground-state bipartite 
entanglement monotone, because for any entangled pure state 
Dschmidi > 2. Indeed, ei coincides with the ground-state bi- 
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partite geometric entanglement, defined as the distance of pc 
from the set of bi-separable pure states \<ps) ® \xr) II 1 211 . On 
the other hand, the fact that ef* is not necessarily a full en- 
tanglement monotone when d > 1 reveals the subtlety of the 
relation between entanglement and frustration. States in a de- 
generate global ground space are entangled just because of the 
linearity of quantum mechanics, and yet this a priori entangle- 
ment needs not be a source of frustration. An exquisite exam- 
ple is the ferromagnetic spin-1/2 Ising model on a triangle. 
This simple model is obviously frustration-free, either clas- 
sically or quantum mechanically (all local interactions com- 
mute). Nevertheless, in the quantum regime, we have that 
any state of the form \tfi) - a\ TTT> + P\ Hi) is an accept- 
able ground state that exhibits a nonvanishing entanglement 
for every possible bipartition as long as a,/3 + 0. The global 
ground space is thus doubly degenerate. On the other hand, 
being the local ground space twofold degenerate as well, for 
the global ground state it will always be Dschmidt < d = 2. 
Regardless of the values of the superposition coefficients, this 

(2) 

implies el = and therefore the absence of entanglement- 
induced sources of frustration. 

According to Eq. (0, the ground state can belong to three 
different classes ]9|]. It will be a frustration free (FF) state if 
fs = 6y" = VS . It will be an inequality saturating (INES) 
state if fs = ef* VS ; therefore, an FF state is a particular case 
of an INES state. Finally, it will be a non inequality satu- 
rate state (non-INES) in all the other situations. As long as 
the systems being considered admit only one non degenerate 
ground state, the same classification applies unambiguously 
also to the corresponding models. 

New issues arise if we consider systems possessing differ- 
ent degenerate ground states. In this case the quantification 
of frustration provided by Eq.dQ will yield in general state- 
dependent results within the same class of models and sym- 
metries. For each model frustration becomes a local, state- 
dependent concept, according to the value taken by Eq.© on 
each different degenerate ground state. On one side, this fea- 
ture confirms that entanglement is a necessary but not suffi- 
cient ingredient to characterize and quantify frustration glob- 
ally in the quantum domain. On the other hand, it implies 
that the measure defined in Eq.(Q~|i should not be applied sepa- 
rately to each degenerate pure ground state. Rather, it should 
be evaluated on an appropriate, average ground state so de- 
fined that it contains all the possible information about the 
global ground space of the system. In this way it will be pos- 
sible to quantify the global frustration properties of a model 



Hamiltonian rather than just its local features in the different 
degenerate ground states. 

In the following, we will show that identifying the ele- 
ments, beyond entanglement, that are needed to characterize 
the global aspects of frustration in the quantum domain, also 
determines the requirements for a state-independent, global 
quantification of frustration. To this end, the crucial observa- 
tion is that statistical mixtures of degenerate ground states are 
themselves legitimate ground states. Moreover, when symme- 
tries are conserved, all the different degenerate ground states 
have the same statistical possibility to be realized. Therefore, 
one can introduce as the appropriate global average ground 
state the maximally mixed ground state (MMGS), that is the 
convex combination with equal a priori weights of all the 
possible degenerate ground states. This principle of a pri- 
ori equiprobability guarantees the correct quantification of the 
global frustration properties of quantum Hamiltonians. The 
MMGS acts as a projector on the global ground space and, at 
variance with single pure ground states, it satisfies all symme- 
tries of the Hamiltonian model being considered. We can thus 
classify the different models and their global frustration prop- 
erties with respect to the properties of the MMGS: a model is 
frustration free on average if its MMGS is frustration free; it 
is INES on average if its MMGS is INES; and it is otherwise 
non-INES on average if its MMGS is non-INES. If a system 
admits a unique, non degenerate ground state, then the local 
and on average classifications coincide. In all other cases if a 
model is locally FF (INES) then it is also FF (INES) on aver- 
age, while in general the inverse does not necessarily hold. 

When addressing mixed states, the entanglement-to- 
frustration relation undergoes essential modifications and gen- 
eralizations. The quantity e§ ceases in general to be a bi- 
partite entanglement monotone when computed directly on 
mixed states IU3I1 and must be replaced by its convex roof 

YjP^f (*j#*>tol) • ( 3 ) 
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In Eq. (f3]) the infimum is taken over all the possible convex 
decompositions into pure states li/^XW of the total system's 
ground state: pc = YjkPkWktity^- The quantities E^ R and 
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ei d) always coincide for pure states, and in general differ for 
mixed states. Physically, this difference comes about because 
the noise present in the reduced density matrix of a globally 
mixed state is associated not only to the presence of entan- 
glement, as for pure states, but also to the presence of classi- 
cal statistical correlations that emerge after local generalized 
measurements on one subsystem Ml44l6ll . 

Total correlations (quantum plus classical) in mixed states 
are usually quantified in terms of the von Neumann entropy; 
however, the relation holds in general: to every entanglement 
monotone there corresponds a type of classical or quantum 
correlation emerging when a local generalized measurement is 
implemented 01711 . For our purposes, we need to evaluate the 
classical correlations between two subsystems, say a and b, in 
terms of e a d \ Although the function is not a full mono- 
tone, it is anyway a proper quantifier of local mixedness, and 
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as such detects correlations between subsystems. Consider 
then a composite systems (made of the subsystems a and b) 
in a mixed state p, and its reduced state p a = ti^p. A general- 
ized measurement on b is defined by a set of positive operators 
{M b (x)}, such that ^ x M b {x) = 1 OH- The measurement 
detects the result x with probability p(x) = tr (1 ® M/,(x)[p]), 
leaving the system in the state p(x) = l a ® M/,(x)[p]. Re- 
placing the von Neumann entropy wit h in the general ex- 
pression for classical correlations IU4[|l5ll . one has: 
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min V p{x)e a d) (tx b p(x)) . (4) 
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The classical correlations C^j are expressed as the difference 

between the total correlations e a d \ evaluated on the reduced 
state p a , and the smallest convex combination of total corre- 
lations, obtained as the minimum over all possible local gen- 
eralized measurements \Mh(x)} on subsystem b. On the other 
hand, it is well known that any n-partite mixed quantum state 
p can be obtained from an (n + l)-partite pure state by trac- 
ing out the degrees of freedom of an ancillary party A. This 
fact allows to generalize to the case of the geometric quan- 
tity el the results originally obtained for the von Neumann 
entropy II 1711 . according to the following 
[Theorem 1 - Purification, entanglement, and classical cor- 
relations]: Given a pure tripartite state \iffsRA) s < Ks® 1 Hr®'Ha 
and its reduced density matrices p SR = ^a\^sra){^sra\, Psa = 
^r\^sra){^sra\ and p s = trgA^SRAi^ssAV then: 
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that is, the total correlations expressed by the geometric quan- 
tity ef 1 are the sum of the bipartite entanglement (convex 
roof) between S and R and the classical correlations between 
S and A. The proof of the theorem is provided in the supple- 
mentary material 02OI1 . We will apply Theorem 1 to the purifi- 
cation of the MMGS with the ancilla A, where the presence of 
classical correlations detected by local generalized measure- 
ments on A is a consequence of the degeneracy of the ground 
state. Conceptually, the ancillary party A can be thought of 
as a suitable quantum reservoir entangled with the bipartite 
system (S \R), yielding the MMGS as the on average reduced 
equilibrium state [21]. Equivalently, A can be seen as a quan- 
tum reference system, such that if one traces over A, every 
ground state of {S \R) is equiprobable, as identified by a com- 
plete set of superselection rules Il22ll . 

Comparing Eqs. (0 and (f2]i we obtain a unified lower bound 
to frustration encompassing both the non degenerate and de- 
generate cases: 
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The above exact inequality yields that, in the quantum do- 
main, frustration is not only due to the underlying geome- 
try, as for the classical case, and/or to entanglement, as for 
systems with non degenerate ground states. In general, it 
depends on the interplay of these two elements with a third 



source, namely statistical correlations established outside of 
the (i-fold degenerate local ground space due to the degener- 
acy encoded in the MMGS. When the ground state is pure and 
non degenerate, the general bound Eq. d6} evaluated on the 
MMGS reduces to Eq. (f2j), as the classical correlations C^ A 



vanish and the convex roof E^ R coincides with e' s . 

Let us illustrate with a simple but nontrivial example the 
difference between the local and on average characteriza- 
tions of frustration provided, respectively, by Eq. (f2) and 
Eq. (O. Consider a ring of five spins with periodic boundary 
conditions described by a ferromagnetic Ising Hamiltonian 
H = - Zf=i S-Sf +1 (S% = S°Va= x,y,z). All local terms in 
the Hamiltonian commute: in this sense, the model is clas- 
sical. The global and local ground states are both two-fold 
degenerate (spin-flip symmetry). Since every element of the 
global ground space is FF, the model is locally always FF, i.e. 
FF on each of the different degenerate global ground states, 
and hence FF on average, that is FF on the MMGS. Next, 
let us modify the Hamiltonian from classical Ising to quan- 
tum XY: H=- Zti(S-Sf +l + AS)S] +l ) , while the geome- 
try of the system remains unchanged. One would now ex- 
pect that frustration in the system should arise from the non- 
commutativity of the local terms in the XY Hamiltonian (now 
the local ground space degeneracy d = 1). Accordingly, one 
verifies that the model is INES on average, i.e. the inequal- 
ity Eq. © is saturated by the MMGS, with the actual values 
of the frustration measure fs and of the total correlations e£ 
depending on the anisotropy A. For instance, for A = 0.1, 
one has fs = e^ 1 ' 0.476 for each of the five different spin 
pair-interaction terms. However, regardless of the value of A, 
such model is never locally INES. Indeed, given the doubly 
degenerate global ground space, let us pick e.g. the ground 
state that is eigenstate of the parity operator along the x di- 
rection with eigenvalue +1. For this pure ground state, the 
measure of frustrations fs takes always the same values as 
for the MMGS (this is actually an interesting general prop- 
erty for all twofold degenerate ground states [20]). However, 
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the bipartite geometric entanglement E S ^ R is always well be 
e^. 1 '. For instance, with A = 0.1 



one has E^ R * 0.001 <sc 4" * 0.476. The impossibility to 



saturate inequality Eq. (f2]i is related to the fact that the se- 
lected pure ground state breaks the symmetry of the model 
Hamiltonian. Modifying the geometry, e.g. by adding a direct 
antiferromagnetic interaction between the first and the third 
spin: H' = S\Si + ASjSo, introduces a further, geometric 
source of frustration. In this case the MMGS no longer satu- 
rates inequality Eq. © and the system ceases to be INES, thus 
signaling the presence of geometric frustration. 

From the above discussions it follows that it would be de- 
sirable to identify a set of conditions to detect a priori the 
frustration properties of the global ground space of a given 
model or class of models. These conditions should include as 
a particular case the ones previously determined for non de- 
generate ground states [9]. Extending the Toulouse criteria {3] 
to the quantum domain, we need to identify a prototype model 
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that is INES on average and then define a group of local oper- 
ations under which the property of being INES on average is 
preserved. We define the prototype model as 
[Prototype model]: A quantum spin Hamiltonian of the type: 



(7) 



is a prototype model if there exists at least one local ground 
space common to all pair local interactions hij and every lo- 
cal coupling vector /y has non-negative components. Having 
defined the prototype model, we state the following: 
[Conjecture 1 - INES property and prototype models]: All 
prototype models are INES on average. 
[Conjecture 2 - INES property and local transformations]: 
Every model obtained from a prototype model by local unitary 
operations on each spin and partial transposition on any ar- 
bitrary set of sites {K} is still INES on average. 

It is evidently quite hard to prove these two conjectures in 
all generality. In the supplementary material [|20|, we pro- 
vide an analytical proof for the one-dimensional quantum 
XY model in the thermodynamic limit, and strong numeri- 
cal evidence obtained for more than 2 • 10 s randomly gen- 
erated models with exchange interactions on arbitrary random 
graphs with a total number of sites N < 9. Preservation of 
the INES property on average under partial transposition is 
to be expected because this property is directly related to the 
absence of geometric frustration in the model. Nevertheless, 
the preservation of the INES property on average is far from 
trivial. Indeed, contrary to what happens under local unitary 
transformations, the properties of degenerate ground states 
of a given Hamiltonian before and after partial transposition 
need not be related. As a straightforward example, let us con- 
sider a spin- 1/2 Heisenberg chain (open boundary conditions) 
of N = 4 spins with homogeneous nearest-neighbor ferromag- 
netic couplings. The global ground state is five-fold degener- 
ate while the local ground space is three-fold degenerate and 
for any couple of interacting spins we have f$ = ef^ = 0. 
Performing partial transposition on spins 2 and 4 (or, equiva- 
lently on spins 1 and 3) the original model maps in a antiferro- 
magnetic Heisenberg chain possessing both a nondegenerate 



global and local ground state with f$ 



0.067 for both 



the pairs of spins (1,2) and (3,4), while f$ = e, 



(i) 



1/2 for 



the pair of central spins (2,3). The transformed model has 
very different ground-state properties compared to the initial 
one, and yet it remains INES. 

In conclusion, we have derived an exact lower bound on 
a universal measure of frustration in the general case of de- 
generate ground states. The bound is expressed as the sum of 
two contributions, one due to bipartite ground-state entangle- 
ment and one due to bipartite classical correlations that are es- 
tablished after local generalized measurements on a quantum 
reservoir or a quantum reference frame. This further source 
of frustration adds to geometry and entanglement, yielding 
a rather complex structure that involves several fundamental 
concepts of quantum physics: entanglement of mixed ground 



states, classical correlations arising from quantum degener- 
acy, and the purification of ground states via quantum reser- 
voirs or quantum reference frames. 

We have showed that the frustration properties of quan- 
tum many-body Hamiltonians are encoded in the maximally 
mixed ground state (MMGS), that is the convex combina- 
tion with equal coefficients of all the degenerate pure ground 
states. Given such on average, global classification, we have 
determined the sufficient conditions for a quantum spin sys- 
tem to achieve the bound, generalizing the results obtained in 
the case of models with non degenerate ground states H. For 
systems with doubly degenerate ground states we have proved 
rigorously that local and average frustration coincide: all de- 
generate ground states, and therefore also the MMGS, exhibit 
the same frustration properties [20]. 

The fact that the residual classical correlations after local 
generalized measurements are identified as a novel source of 
frustration in quantum many -body systems may open interest- 
ing research insights concerning competitions among quantal 
and non-quantal aspects in Hamiltonian spin dynamics, as in 
the case of the anisotropic Heisenberg models in an external 
field, or competition between quantal dynamics and thermal 
fluctuations. Future investigations should address the role that 
the presence/absence of the INES property on average actu- 
ally plays in computational 112311 , information- theoretic 
12711 . and thermodynamic characterizations 11281 12911 of many 
body quantum systems. Experimental consequences might 
soon be derived and tested, as the controlled quantum simu- 
lation of magnetism, classical and quantum, already includes 
the first demonstration of antiferromagnetic spin chains with 
optical lattices ll30Tl . the probing of small frustrated Heisen- 
berg spin systems by photonic simulators 03111 . the realization 
of classically frustrated Ising spins with trapped ions 113211 and 
optical lattices 13311 . up to the recent comprehensive proposals 
for the quantum simulation of large classes of quantum frus- 
trated magnetism with ion crystals [34] and color centers in 
diamond 113511 that might open the way to the precise verifica- 
tion of long-standing predictions on exotic phases of matter. 
In this respect, the use of tools inspired by quantum informa- 
tion science can lead to a broader and deeper understanding 
of collective quantum phenomena via long-range "entangle- 
ment patterns" ir36Tl . beyond the Landau-Ginzburg framework 
of symmetry breaking and local order parameters [37 1. 

F.I. and S.M.G. acknowledge financial support through the 
FP7 STREP Project iQIT, Grant Agreement No. 270843. 



SUPPLEMENTARY MATERIAL 



MONOGAMY OF CORRELATIONS 

We begin by reviewing the consequences of the Cauchy in- 
terlacing theorem [ 10] with respect to convexity and concavity 
of entanglement monotones, either faithful or unfaithful. 
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Proposition 1 (Concavity of the generalized geometric en- 
tanglement). The generalized geometric entanglement e^ 1 is 
concave: 



ef(p)>Y J p k ef(p k ), 



(8) 



provided p = ZkPkPk ■ 



Proof. From Cauchy interlacing theorem we have that, for all 
k and for all orthonormal bases 

d d 

If \\i))i=i,...,d are the eigenvectors corresponding to the first d 
largest eigenvalues {/I- }i=i,...,d of p, ordered in descending or- 
der, one has that: 

d d 
i=l i=l 

d 

= £j2^Pk(i\Pk\i) 

k i=l 
d 

< Yu Z M a Cp*) • 



p x = tr[(lj ®M x )psa\, leaving the total system in the state 



Px = j- x ds ®M x )p SA . 



A bona fide measure of the classical correlations C ( ^ A exist- 
ing between the two subsystems of globally bipartite system, 
satisfying a minimal set of axioms, was defined as the dif- 
ference between the total correlations in the global state p$A, 
measured by the von Neumann entropy of the reduced state 
ps and the residual quantum correlations after a local gener- 
alized measurement on party A B14. Ilal . Here we introduce 
an analogous definition, except for the crucial difference that 

the von Neumann entropy is replaced by the concave function 

id). 



c (d) = ■ 
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(p s ) - min V p x ef (tr A p x ) . (11) 



k i=l 



As we will now show, this is still a bona fide measure of 
bipartite classical correlations that satisfies the Koashi-Winter 
monogamy relation originally proven for the entanglement of 
formation II 1711 : 

Proposition 2 (Koashi-Winter monogamy of correlations). 

Given a pure three-partite state \iffsRA) e "Ms ® 'Mr <S> 'Ha 
and their reduced density matrices psr = ^a\i//sra)(^sra\, 
Psa = ^r\^sra){^sra\ and p s = tr ra\<J'sra)(<{'sra\, the fol- 
lowing relation holds 



(d) 



On the other hand, every bipartite mixed state psr of a bi- 
partite system (S \R) can be decomposed into infinite convex 
combinations or pure states: 



PSR = Y_ l PkWk){fk\, Pk>0, = 1 



(9) 



and a faithful measure of the entanglement in p$R can be de- 
fined via the convex roof construction over pure-state entan- 
glement: 



k 



(10) 



where the infimum is taken over all the possible convex de- 
compositions into pure states, Eq. (TT0T >. 

The degree of mixedness in the reduced density matrix p$ 
of a bipartite mixed state psr is directly related not only to 
the degree of bipartite entanglement of psr but also to a set 
of classical statistical correlations, namely all the correlations 
that are left after a local generalized measurement is per- 
formed on one party (subsystem) of the total system lfl4[|l5ll . 
Given a bipartite system (S\A), a generalized measurement 
(POVM) on party A is described by a set of positive operators 
{M x }, such that £ A . M x = 1 HI 8IT . In a bipartite mixed state psa, 
with subsystem S in the reduced state ps - tr^p^, a general- 
ized measurement on A yields some result x with probability 



Md) 
'S\A 



(12) 



Proof. An analogous proposition was originally proved 
in 01711 . with the entanglement of formation replacing the ge- 
ometric monotone of entanglement E^ R defined by Eq. ( [Tol l, 
and the standard classical correlations introduced in Refs. [ 14, 
[l5ll replacing the classical correlations defined by 

Eq. (fTTT >. All these differences reduce to the function 
replacing the von Neumann entropy, used throughout in 
Ref. 11*711 . On the other hand, the only property of the von 
Neumann entropy which enters in the proof for the original 
proposition of Ref. lfl7ll is concavity. Since the function 
is concave, as proved in Proposition [I] it follows that Propo- 
sition|2]holds. □ 

If the state psr is pure, then the classical correlations C^ A 
vanish, and the total correlations reduce to the generalized 
geometric entanglement Eg) R . If psr is mixed, there are many 
possible purifications of which psr is the reduced state 

upon tracing out the degrees of freedom of party A. This fact 
implies in particular that: 

Proposition 3. The classical correlations C ( ^ A do not depend 
on the purification procedure. Equivalently, all purifications 
\i//sra) coincide up to a local unitary transformation on party 
A. 

Proof. In order to prove this property explicitly, one needs 
first to recall that, by construction, the classical correlations 
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C S ' A , just like entanglement, are invariant under local unitary 
operations [14]. Moreover, let us recall that purification of 
a mixed state p SR amounts to attaching an orthonormal basis 
{\af)}j to its eigenvectors {\v ( f R )}j as follows: 

j 

i 

If {Ifl^')}; were not an orthonormal basis, the reduced state 
would not coincide with psr. If psu does not possess degen- 
erate eigenspaces, then the only freedom is on the choice of 
the basis {\af)}j, which is equivalent to performing a local 
unitary operation on party A. If psr possesses some degen- 
erate eigenspace, there is freedom in the choice of its eigen- 
vectors {\v ( f R }}j. One can purify the degenerate eigenspace 
spanflv™) \jtP] of psr with eigenvalues p, in the following 
ways: 

jeP jeP 
jeP 

l*£>=Z&»|v&>, d4) 

leP 

where %a are the matrix element of a unitary matrix, in order 
for to be orthogonal vectors. Thus 

jeP leP 

and the freedom of purification reduces again to a unitary op- 
eration on the party A. This local unitary freedom can be con- 
sidered as a part of the POVM {M x }. Since the minimization 
over all the POVMs enters in the definition Eq. (fl~T), the clas- 
sical correlations C^ A are independent of the purification. □ 

State-independent frustration in systems with doubly 
degenerate global ground states 

Consider spin- 1/2 models described by Hamiltonians of the 
type 

H =~Y. J ij S i S r < 16 > 

where the operators are the usual spin operators along the 
fi direction acting on the i-th spin. Such Hamiltonians are in- 
variant under the action of the parity operator P a = Ylj °" a 



(with a - x,y, z), that is [H, P a ] = 0. On the other hand, 
the parity operators along different directions do not commute 
when the system is made of an odd number of spins. This fact 
implies that the global ground state of H cannot have degener- 
acy dc < 2. In other words, the global ground state cannot be 
unique. In the class of Hamiltonians Eq. ( fTo*b we consider the 
case in which the global ground state is precisely two-fold de- 
generate. This instance is realized by the XY model, either in 
a finite chain with an odd number of spins l38ll or in the ther- 
modynamic limit ll39tl . We emphasize that higher degenera- 
cies (dc > 2) may occur either for models such that there are 
strong resonances in the Hamiltonian parameters, and/or such 
that they obey additional symmetries, so that indeed twofold 
degeneracy (dc = 2) holds typically. According to the afore- 
mentioned invariance under parity transformations, we can in- 
troduce as a basis for the global ground space the two eigen- 
states of the parity operator along a fixed direction a with 
eigenvalues ±1. These two eigenstates are denoted, respec- 
tively, as the even and odd states and |— a ). Up to a global 
phase, one has that |+ ff ) = Pp\± a ), with ft _L a. As a conse- 
quence, the measure of frustration f$ coincides on both states 
for every spin pair. All pure global ground states can be writ- 
ten as = a\+ a ) + b\- a ), with \a\ 2 + \b\ 2 = 1. Given a spin 
pair S and n,s the projector on any Bell state defined on S , 
one has that 

tr(|^|n s )=tr[(M 2 |+«X+ ff | + \b\ 2 \- a )(- a \)n s ] 

+a*bxi(\+ a )(- a \n s ) + c.c. . (17) 

Recalling the invariance of the Hamiltonian under the parity 
operations and that |— „) = Pp\+ a ), one has further that 

tr(l+„X+*|n s ) = tr(K><- ff |n s ) . (18) 

To evaluate terms in the r.h.s. of Eq. [l7]one can re-write the 
states \± a ) in the relative state representation [40] with respect 
to the Bell basis of spin pair S . Taking into account that every 
Bell state is an eigenstate of some parity operator, one has that 

\± a )=Yj (gWrehWSh+Glip.hVBSh) . (19) 

£=±1 

where s runs on the possible eigenvalues {± 1 } of the parity op- 
erators P a acting on subsystem S , \i// e )s and \<p e )s are the two 
Bell states of spin pair S such that when P a is applied on them 
it yields the eigenvalue e, while \^^)r and \^)r are states de- 
fined on the remainder R of the entire system. Since |± tr ) are 
eigenstates of P a , also I^Xr and \0^)r must be eigenstates of 
the same parity operator when acting on the Hilbert space of 
R. Moreover, because \+ a ) and yield two different values 
of the parity, we must have that |Tg )r (or IO^Xr) must yield a 
different value of the parity than \^ )r (or \<&~)r). Collecting 
these simple facts it is straightforward to verify that 

trains) = tr(|+„X+„|ns) . (20) 

Recalling the definition of the measure of frustration f$, 
Eq. (1) in the main text of the present letter 114 ill , from Eq. (|20| ) 
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it follows that the local frustration, that is fs evaluated in any 
degenerate pure ground state, does not depend on the chosen 
ground state in the case of a doubly degenerate ground space. 
Moreover, since the measure f s is linear, the global, on aver- 
age frustration of the maximally mixed ground state (MMGS), 
that is the convex combination with equal weights of all the 
pure degenerate ground states, is the convex combination with 
equal weights of the local frustrations of the pure degenerate 
ground states. Since we have just seen that for doubly degen- 
erate ground spaces the local measures of frustration coincide 
for all pure degenerate ground states, it follows that the global 
and local frustrations coincide. 



EVIDENCE SUPPORTING CONJECTURES 1 - INES AND 
PROTOTYPE MODELS; AND CONJECTURE 2 - INES AND 
PARTIAL TRANSPOSITION 

Numerical results for random graphs with a finite number of 
spins 

Here we will discuss in detail the verification of the two 
Conjectures put forward in the main text of the present let- 
ter 1141 1] , that is Conjecture 1: All prototype models are INES 
on average, and Conjecture 2: Every model H' = H Tk ob- 
tained from a prototype model H by partial transposition on 
any arbitrary set of sites {K} is still INES on average. 

In order to gain as much insight as possible on the validity 
of these two Conjectures, we have analyzed large numbers of 
randomly generated Hamiltonians without geometric frustra- 
tion, and checked that they satisfy the equality fs = ei for 
all pairs of spins (subsystems) S entering in the local inter- 
action terms h$ . Instances of the prototype model have been 
generated, where the gauge was fixed at the unitary sector of 
the gauge group, since this symmetry is trivially proven an- 
alytically. In order to check Conjecture 2, a random gauge 
transformation consisting of pure parity transformations was 
applied. In total 2 x 10 5 geometric frustration-free Hamilto- 
nians have been checked by exact diagonalization. No single 
instance was found that fails to satisfy fs = ei® for each pair 
of spins. 

The numerical tests were performed according to the fol- 
lowing route. We have generated random models of the XYZ 
(general anisotropic gapped spin- 1/2 models with exchange 
interactions), XXZ (gapless Heisenberg) and XXX (gapless 
isotropic Heisenberg), satisfying the conditions described in 
the main text 14111 . either with homogeneous or inhomoge- 
neous nearest-neighbor couplings. The models have been gen- 
erated first by generating a connectivity graph (fV, £) with 
\ r V\ = n vertices connected by edges <5. A spin-1/2 opera- 
tor was placed at each vertex. For each edge (if) e S, an 
interaction of the XYZ, or XXZ, or XXX type was generated. 
Depending on whether an homogeneous or inhomogeneous 
model was produced, the interactions were either generated 
by replicating them identically for all edges, or all the single 
interactions were generated independently. Finally, either the 



identity or a parity transformation were randomly applied to 
every spin, each with probability 1 /2. After having generated 
randomly both the lattice and the model Hamiltonian, we have 
diagonalized numerically the different Hamiltonians and for 
all pairs of nearest-neighbor interacting spins we have deter- 
mined the reduced density matrices and evaluated the measure 
of frustration fs and th geometric quantity e s d \ All the numer- 
ical tests have been performed programming with MATHE- 
MATICA 7. 

In Tables U HI] and [Till we summarize the results of the 
extended numerical analysis, reporting the total number of 
models generated for each type of interaction, the number of 
Accepted Hamiltonians (those for which the INES condition 
fs = ei is satisfied for all bonds (if) e £), and the number of 
Rejected Hamiltonians (those for which the INES condition 
fs = is violated at least for some bond). The total num- 
ber of spins (vertices) ranged from 5 to 9. The conclusion is 
that no rejected Hamiltonian has been found, meaning that no 
violation of Conjectures 1 and 2 could be identified. 

n Accepted Rejected 
~ 30000 
7 20000 
9 10000 

TABLE I: Results for XYZ models. 

n Accepted Rejected 

5 30000 
7 20000 
9 10000 

TABLE II: Results for XXZ models. 

n Accepted Rejected 

5 30000 
7 20000 
9 10000 

TABLE III: Results for XXX models. 



Analytical results for the XY chain 

Our conjectures can be proved analytically for the case of 
the XY quantum spin chain in the thermodynamic limit, ex- 
ploiting the exact solution available for the model [42], whose 
Hamiltonian reads 

H *r= \ Tl 1 + y)S ^ + (1 " ' (21) 

i 

where y is the anisotropy parameter that can vary in the inter- 
val [0,1], realizing the standard (classical) Ising model (y = 
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FIG. 1: Behavior of the bipartite geometric entanglement (black 
solid line) and of the universal measure of frustration f s (red dots) 
as functions of the anisotropy parameter y for a pair of neighboring 
sites in the XY quantum spin- 1/2 chain in the thermodynamic limit. 



1), the generic (quantum and gapped) XY model (0 < y < 1), 
and the gapless, critical XX model (y = 0). Except for the 
Ising case (y = 1) the local interactions between neighbor- 
ing spins allow to accommodate a unique non degenerate lo- 
cal ground state while the global ground state is doubly de- 
generate 11381 13911 except for the XX point (y = 0) at which 
the global degeneracy is removed. Throughout, we consider 
the thermodynamic limit of XY chains with periodic bound- 
ary conditions, so that translational invariance is enforced at 
all stages. In Fig. Q] we thus report the behavior of the bipar- 
tite geometric entanglement ei as a function of y (black solid 
line) for an arbitrarily chosen pair S of neighboring spins and 
we compare it with the measure of frustration fs (red dots) for 
< y < I. One sees that the two quantities coincide exactly 
for all values of the anisotropy parameter, proving exactly our 
two Conjectures for this class of models. In the plot, the point 
y = 1 is excluded. In fact, at this point the local ground state 
becomes doubly degenerate and one must then consider the 



geometric quantity rather than e)!'. At y — 1 the XY chain 
turns in the classical, one-dimensional, frustration-free Ising 
model on a infinite chain. Indeed, it is straightforward to ver- 
ify that both e s and fs vanish identically for this case. 
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